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Abstract 

We attempt here to understand successfully some crucial aspects of J/^'-production in 
some high energy nuclear collisions in the light of a non-standard framework outlined in the 
text. It is found that the results arrived at with this main working approach here is fairly in 
good agreement with both the measured data and the results obtained on the basis of some 
other models of the 'standard' variety. Impact and implications of this comparative study 
have also been precisely highlighted in the end. 
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1 Introduction 



The study of the J/^'-mesons in ultra-relativistic heavy ion colhsions has consistently been 
considered to be a potentially powerful tool for studying the properties of the hypothetical 'hot 
and dense matter' created in these collisions. Predictions about the suppression of J/^'-meson 
production in the nuclear collision, at certain stages [Ij and an anomalous enhancement at certain 
other stage as well [2] are always treated as very powerful diagnostics. The deciding factor, in 
terms of the Standard Model (SM), in both the cases, is the number of charm-anticharm (cc) 
pairs (Ncc) created in the early stage of hard parton collisions in A + ^ (or in A + B) reactions. 
On the basis of it, one might arrive at the assumptions of the J /^-suppression (at Ncc < 1) and 
the J /^-enhancement (at Ncc > 1) effects. 

In the past we had constantly refrained from giving such undue importance to any 

controversy about either suppression/enhanced suppression or an enhancement of J/^ mesons 
in reality. Rather our primary intent would centre around understanding and interpreting the 
nature of data of inclusive cross-sections and some other important observables of J/^ mesons in 
BNL-RHIC and CERN-LHC experiments. In fact, we have an altogether different and heretic 
view (which will be outlined very briefly in the next section) about the mechanism of J/^ 
production in high energy particle and nuclear collisions. 

We define our objectives here as: 1) to explain the main and major features of the latest data 
on J/^-production in BNL-RHIC and CERN-LHC experiments from the proposed alternative 
approach built up by us in a set of previous works done in the both the remote and recent past 
[6]-[7] and 2) to compare our model-based calculations with some other competing models. 

Our plan of work presented here is as follows. In section 2, we provide brief outlines of the 
models chosen for this study. The section 3 provides very brief outlines of the models which are 
founded on the gantlets of assumptions of the 'Standard'-model and which have been reckoned 
with for the sake of comparison. In section 4 we give the results and general discussion. And in 
the last section (section 5) we offer the final remarks-cum-conclusions. 

2 The Main Approach: An OutUne 

In the present work, we will make use of a theoretical model to interpret some of the latest and 
topical observables of the J/^* production which were measured and reported by the different 
groups in the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) exper- 
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iments in the recent past. The model has a modest degree of dynamical basis and some prior 
check-ups with data [3]-[5]. It is called here as the sequential chain model (SCM). 

The outline and features of the model we use here are obtained, in the main, from our previous 
works [3j-[5]. According to this model, high energy hadronic interactions boil down, essentially, 
to the pion-pion interactions; as the protons are conceived in this model as p=(tt^tt^i!}), where i? 
is a spectator particle needed for the dynamical generation of quantum numbers of the nucleons. 
The multiple production of J/^'-mesons in a high energy proton-proton collisions is described 
in the following way. The secondary vr-meson or the exchanged ^^-meson emit a free w-meson 
and pi-meson; the pions so produced at high energies could liberate another pair of free g and 
trapped w-mesons (in the multiple production chain). These so-called free g and w-mesons decay 
quite a fast into photons and these photons decay into ^ or ^' particles, which, according to 
this alternative approach is a bound state of ^lO. or QTI' particles. 

The inclusive cross-section of the ^'-meson produced in the pp collisions given by 



d^a 1 -5.35(p^ + m],^) 

^TT\p+p^J/'f+x - Cj/*Tv7^^P(^ ^r~n r)exp(-1.923< nj/^ >ppx), (1) 

dp pj.^ < nj/^j, >pp (I - X) 



where the expression for for average multiplicity for \I'-particles in pp scattering would be given 

by 

< nj/* >pp = 4 X lO-S^/^ (2) 

In the above expression, the term |Cj/,j,| is a normalisation parameter and is assumed here to 
have a value = 0.09 for Intersecting Storage Ring(ISR) energy, and it is different for different 
energy and for various collisions. The terms px, x and mj/^ represent the transverse momentum, 
Feynman Scaling variable and the rest mass of the J/^ particle respectively. Moreover, by 
definition, x = 2pl/^/s where is the longitudinal momentum of the particle. The s in 
equation (2) is the square of the cm. energy. 

The second term in the right hand side of the equation (1), the constituent rearrangement 
term arises out of the partonic rearrangements inside the proton. It is established that hadrons 
(baryons and mesons) are composed of few partons. These rearrangements mean undesirable 
loss of energy , in so far as the production mechanism is concerned. The choice of N^i would 
depend on the following factors: (i) the specificities of the interacting projectile and target, (ii) 
the particularities of the secondaries emitted from a specific hadronic or nuclear interaction and 
(iii) the magnitudes of the momentum transfers and of a phase factor (with a maximum value 
of unity) in the rearrangement process in any collision. The parametrisation is to be done for 
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two physical points, viz., the amount of momentum transfer and the contributions from a phase 
factor arising out of the rearrangement of the constituent partons. Collecting and combining all 
these, we propose the relation to be given by |8| 

Nr = 4.< Npart >'/^ e, (3) 

where < Npart > denotes the average number of participating nucleons and 9 values are to be 
obtained phenomenologically from the fits to the data-points [9]. 

In order to study a nuclear interaction of the type A + B — t- Q + x, where A and B are 
projectile and target nucleus respectively, and Q is the detected particle which, in the present 
case, would be J /"if -mesons, the SCM has been adapted, on the basis of the suggested Wong fW\ 
work to the Glauber techniques. The inclusive cross-sections for J/^ production in different 
nuclear interactions of the types A + B ^ J/^ + X in the light of this modified Sequential 
Chain Model (SCM) can then be written in the following generalised form as: 

E-^\a+b^j/^+x = aj/^PT~'^'^ exp(-c(p^ + "ij/vE,))exp(-1.923< rij/^ >pp x). (4) 

where aj/^, Nji and c are the factors to be calculated under certain physical constraints. The 
set of relations to be used for evaluating the parameters aj/^ is given below. 

_ 3 {AaB + BaA) 1 

Here, in the above set of equations, the third factor gives a measure of the number of wounded 
nucleons i.e. of the probable number of participants, wherein AaB gives the probability cross- 
section of collision with '5' nucleus (target), had all the nucleons of A suffered collisions with 
-B-target. And BaA has just the same physical meaning, with A and B replaced. 

Besides, in expression (5), the fourth term is a physical factor related with energy degradation 
of the secondaries due to multiple collision effects. The parameter a' occurring in eqn.(5) above 
is a measure of the fraction of the nucleons that suffer energy loss. The maximum value of a' is 
unity, while all the nucleons suffer energy loss. This a' parameter is usually to be chosen |10j . 
depending on the centrality of the collisions and the nature of the secondaries. 

The "a" factor in the expression (5) accommodates a wide range of variation because of 
the existence of the large differences in the in the normalizations of the J/^ cross-sections for 
different types of interactions. 
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3 Some Competing Models: The Brief Outlines 

3.1 The Gluon Saturation Approach 

In this approach [llj, it is assumed that the nuclear wave functions in very high-energy nuclear 
collisions can be described by the Color Glass Condensate (CGC). The primary effect is the 
suppression of J/^' production and narrowing of the rapidity distribution due to saturation of 
the gluon fields in heavy ion collisions relative to p + p collisions. In addition, the production 
mechanism is modified from p + p such that the multigluon exchange diagrams are enhanced 
in heavy ion reactions. It should be noted that this model does not include any hot medium 
effects, but does have a free parameter for the overall normalization factor. 

3.2 The Quark Coalescence Model (QCM) 

The Quark Coalescence Model (QCM) p^, is a two-stage simulation. In stage I the incoming 
target and projectile nucleon interactions are tracked, a fluid of pre-hadrons are formed and in 
stage II the produced pre-hadrons interact and decay in a standard relativistic cascade model. 
The time history of all the collisions recorded in stage I sets up the geometry and initial conditions 
for stage II. Basic inputs for the simulation are measured hadron-hadron crosssections, rapidity, 
transverse momentum and particle multiplicity distributions, of which the last are conforming 
to KNO scaling. The basic assumption of the model is coalescence of charm quarks and anti- 
quarks into charmonia (cc — )• charmonium), which can also be viewed as heavy pre-hadrons in 
stage I. Only a small percentage of charm quarks are expected to coalesce into bound charmonia 
i.e. J/^, the remainder of such quarks appear ultimately as open charm mesons. At RHIC 
energies comover suppression of directly produced charmonia is expected to be large, due to the 
increased particle numbers and densities. 

3.3 Double Color Filter-Oriented Approach 

A mechanism called double color filtering for cc dipoles [13j . makes nuclei significantly more 
transparent in AA than pA collisions. The assumption of this model is as follows: In the cm. of 
nuclear collision the nuclear disks passing through each other leave behind a cloud of radiated 
gluons creating a dense matter, which the J/^ propagates through. In this reference frame the 
J/^ full momentum is px, which ranges from zero to several GeV in RHIC data. Such a low 
energy cc dipole develops the J/^ wave function pretty fast, during time tj < 0.5 fm, which is 
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about the time scale of the medium creation. Thus, what is propagating through the medium is 
not a small cc dipole, but a fully formed J/'^ . The observed nuclear effects in J/'^ production in 
AA collisions is interpreted as a combination of final state interaction (FSI) of J /'^ in the dense 
medium, and the initial state interaction (ISI) effects in production of J/'^ caused by multiple 
interactions of the colliding nuclei. 



4 The Results 

Now let us proceed to apply the chosen model to interpret some recent experimental results of 
J/^'-production reported by various groups for different collisions like p + d + Au, Cu + Cu, 
Au+Au at RHIC and p+p, Pb+Pb at LHC. Here, the main observables are the invariant yields, 
rapidity distributions and the nuclear modification factors which would come under purview of 
the present work. 

4.1 J/\& (total) Crosssections and Rapidity Distribution in p + p Interactions 

As the psi-productions are generically treated rightly as the resonance particles, the standard 
practice is to express the measured J/^ (total) crosssections times branching ratio to muon or 
electrons, i.e. for lepton pairs , that is by Bu/ap!^. 

By using expression (4) we arrive at the expressions for the differential cross-sections for the 
production of J/^'-mesons in the mid and forward-rapidities (i.e. \y\ < 0.35 and 1.2 < \y\ < 2.2 
respectively) in p + p collisions at ^/snn=200 GeV at RHIC. 

^^/Z'^^lp+P^j/M = e.lpr'-''' exp[-0.13(p| + 9.61)] for \y\ < 0.35, (6) 

and 

I ^ d^(T . „ . -1.183 



= 6.5p^^-^«^exp[-0.16(p^ + 9.61)] for 1.2 < \y\ < 2.2. (7) 

For deriving the expressions (6) and (7) we have used the relation x ~ ^^^"^ = ^"^t smhyem j24j 
where niT , Vcm are the transverse mass of the produced particles and the rapidity distributions. 
mj/^ — 3096.9 it O.OllMeF [l^ and Bur, the branching ratio is for muons or electrons i.e. its 
for lepton pairs J/^' — )• /i^/i^/e+e^, is taken as 5.93 it 0.10 x 10^^ in calculating the above 
equations. 

In a similar fashion, the total inclusive cross-sections for the production of J/^'-mesons in 
different rapidities (i.e.|y| < 0.75 and 2.5 < |y| < respectively) in p + p collisions at y/SNN=7 
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TeV at LHC would be, 
(fa , 



^^lp+P^J/M = 78.54p^3.i35g^p[_0_013(^2 ^g g;^)] |y|<0.75, (8) 



and 

37.52p;^^-i3^ exp[-0.016(p| + 9.61)] for 2.5 < \y\ < 0. (9) 



In Fig. 1(a) and Fig. 1(b), we have drawn the sohd hnes depicting the SCM model-based 
results with the help of above equations (6), (7), (8), (9) against the experimental measurements 
[15], [16] respectively. 

For the calculation of the rapidity distribution from the equation (4) we can make use of a 
standard relation as given below: 

dN r 1 d^N , , , 

dpT (10) 



dy J 2ttpt dprdy 

For p + p collisions, the calculated rapidity distribution equation at RHIC-energy ^/snn = 200 
GeV is 

^Ip+p^j/vp+x = 1-215 X 10"^exp(-0.23sinhycm), (H) 

In Fig. 2 we have plotted the rapidity distributions for J/^'-production in p + p collisions at 
y^sjvTV = 200 GeV. Data in the figure are taken from Ref. [T7] and the line shows the SCM-based 
output. 

In a similar fashion, the SCM-based rapidity distribution equation at LHC-energy ^/sNN = 
7 TeV in p + p collisions has been given hereunder 

dN 

— \p+p^j/q,+x = 8.025 exp(-0.043sinhycm), (12) 

In Figure 3, we have drawn the solid lines depicting the SCM-based results with the help of 
above equation (12) against the experimental measurements }18j . 

4.2 Invariant Yields in d + Au, Cu + Cu and Au + Au Collisions at RHIC-energy 
= 200 GeV 

From the expression (4), we arrive at the invariant yields for the J/^-production m d + Au — t- 
J/^ + X reactions for mid and forward-rapidities. 

^ d'^^ I _ V OK in-7^-0.629, 



2TrpT dpxdy 



\d+Au^j/^+x = 7.25 X 10-^p^^-^^'^exp[-0.13(p^ + 9.61)] for \y\ < 0.35, (13) 



7 



Cu+Cu^j/^+x = 6.40 X exp[-0.13(p^ + 9.61)] for \y\ < 0.35, (15) 



and 

= 4..25xlO-^p:^''-^^^exp[-0.l6{pl+9M)] for 1.2 < \y\ < 2.2. (14) 

In Figure 4, we have drawn the solid lines depicting the SCM-based results with the help of 
above two equations (13) and (14) against the experimental measurements |19j . 

For the case of Cu + Cu most central collisions (0-20%) at RHIC, the SCM-based calculated 
theoretical invariant yields for the rapidities |y| < 0.35 and 1.2 < \y\ < 2.2 are given by the 
following equations respectively; 

1 d^N 
2'KpT dpxdy 

and 

^^^b«+c«->J/*+x = 6.38 X 10-5p/-«29 exp[-0.16(p|. + 9.61)] for 1.2 < \y\ < 2.2. 

(16) 

The experimental results for the invariant yields of J/^ production as a function of transverse 
momenta for Cu + Cu collisions are taken from Ref. |20] at centrality — 20% and are plotted 
in Fig. 5. The solid lines in the figure show the SCM-induced results. 

Similarly, for Au + Au collisions at y^sjv7v=200 GeV at RHIC, the equations of transverse 
momenta spectra for 0-20% centrality regions are given by the undernoted relations. 

1 d^N 
2'KPT dprdy 

and 

1 d^N , 



Au+Au^j/^+x = 1-32 X 10-^p^i-°23 exp[-0.13(p^ + 9.61)] for \y\ < 0.35, (17) 



2^^^^^^^^,A«+An^j/vi.+x = 0.91 X W~''p^'-"''exp[-0.16{p'T + 9.61)] for 1.2 < \y\ < 2.2. 

(18) 

For calculating the values of Nji, in general, we have used the values of < Npart > from [llj, |19] . 
In the Fig. 6, the solid lines are the plots of SCM-based invariant yields vs. px as described by 
equations(17) and (18) at forward and mid-rapidities for Au + Au collisions, while the dotted 
curve in the Fig. shows results of coalescence model [12]. The experimental data points in the 
Fig. 6 for the invariant yields of J/^ production as a function of transverse momenta at centrality 
values — 20% and at the rapidities |y| < 0.35 and 1.2 < \y\ < 2.2 respectively are taken from 
the PHENIX Collaboration 
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4.3 The Nuclear Modification Factor Rab 

There is yet another very important observable called nuclear modification factor (NMF), de- 
noted here by Raa which for the production of J/'f is defined by [20] 



d^Nf/l,/dpTdy 

" < N^oiiib) > d^Nj^^/dpTdy- ^^^^ 
the SCM-based results on NMFs for Cu + Cu and Au + Au collisions for forward rapidities are 
deduced on the basis of Eqn.(6), Eqn.(15), Eqn.(17) and Eqn. (19) and they are given by the 
undernoted relations 

Raa\cu+Cu^j/^+x = ^A2p^f^\ (20) 

and 

Raa\au+Au^j/^+x = 0.36p?^i^ (21) 

herein the value of < Ncoii{b) > to be used is ~ 170.5 it 11 [21J for Cu + Cu collisions and for 
Au + Au collisions it is taken as ~ 955.4 it 93.6 |22| . 

In Fig. 7(a), we plot Raa vs. pt for — 20% central region in Cu + Cu and Au + Au collisions. 
The solid lines in the figure show the SCM-based results against the experimentally measured 
results [19] , [20] . The dotted lines in the fig. represent the double Color Filtering approach [13] . 
And in Fig. 7(b), we plot Raa vs. Npart for — 20% central region in Au + Au collisions at 
y^s/vTV = 200 GeV and for Ph + Ph collisions at ^/snn = 2.76 TeV . The solid lines in the figure 
show the SCM-based results against the experimentally measured results from Refs. [TT] and 
respectively. The dashed lines in the fig. represent the Gluon Saturation approach |llj . 



5 Summary and Outlook 

Let us first concentrate on what we have achieved here: (i) The features related to pT-spectra for 
J/'^ production in some particle-particle and nuclear collisions at various high energies have been 
reproduced quite successfully, (ii) The characteristics of rapidity spectra in pp collisions at TeV 
energies have been brought out somwhat modestly satisfactorily, (iii) The features of nuclear 
modification factors in Cu+Cu, Au+Au and Pb+Pb reactions at RHIC and LHC energies have 
been worked out with the help of the applied model. Besides, some of our model-based results 
have also been compared with the performances on the same observables by some competing 
models of 'standard' variety. And these comparisons with data and the results obtained by 
some other models reveal that SCM-based results are in better agreement with data than the 
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other approaches grounded on the 'Standard' model ilk. And this observation is no accident or 
coincidence. In the past such were the recurrent observations made by us valid for many other 
obscvablcs measured in the various high statistics liigh energy particle and nuclear experiments. 

Thus, summing up our past experiences and considering the weightage of the results reported 
here, we are forced to comment finally that this work essentially represents a case of paradigm 
shift in the domain of particle theory, as we have eschewed the conventional views of cc approach 
to J/* production in the 'standard' framework. And this is just the reflection of our radical 
views about the particle structure and the nature of particle collisions. Obviously we obtain the 
fair agreement with data on some observables without inductions of (i) any QGP concept, (ii) 
any prognosis of suppression or enhancement of J/\l'-production. The production characteristics 
resemble all other hadrons. 
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Figure 1: Plot of the invariant cross-section for J/4' production in proton-proton collisions at 
(a) \/s^7v ~ 200GeV and (b) -v/s^vjv ~ 7TeV as function of px- The data points are from [15] 
for (a) and from [TB] for (b) . The sohd curves show the SCM-based results. 




Figure 2: Plot of the rapidity distribution for J/^ production in proton-proton collisions at 
V^NN ~ 200Gey as function of y. The data points are from [19], [T7]. The solid curves show 
the SCM-based results. 
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Fi gure 3: Plot of the rapidity distribution for J /"^ production in p+p collisions at ^/s j^j^ — 7TeV 
as function of y. The data points are from [18j. The solid curves show the SCM-based results. 
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Fi gure 4: Plot of the invariant yields for J /"^ production in d-\-Au collisions at \/sjva'" — 

2mGeV 

as function of px- The data points are from [19j. The solid curves show the SCM-based results. 



o 



Cu+Cu=J/1'+X 
0-20%x5(|y|<0.35) 
SCM 

0%x5(1.1<|y|<2.2) 




2 4 6 

Py(GeV/c) 



Figure 5: Transverse momenta spectra at \y\ < 0.35 and 1.2 < \y\ < 2.2 for J/^ production in 
Cu + Cu central collisions at -v/s^tv ~ 200GeV. The data are taken from [20]. The solid curves 
depict the SCM-based results. 
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Figure 6: Plot of the invariant yields for J/^ production in Au + Au collisions at -v/sjvTV ~ 
200GeV as function of pT- The data points are from ^17j. The solid curve shows the SCM-based 
results while the dotted curve depicts the Coalescence Model |12j . 
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Fi gure 7: (a)Plot of the Raa for J/^ production in Cu-\-Cu anAAuA-Au collisions at \fsj^j^ — 
200GeV as function of px- The data points are from [19], [20]. The solid and dotted curves 
show respectively the SCM and the Double Color Filter-oriented [13] results, (b) Plot of the 
Raa vs. Npart for J/^ production in Au + Au collisions at -v/s^^ = 200GeV and Pb + Pb 
collisions at = 2.76TeV. The red circles depict Au + Au collisions [11 j while the black 

squares represent Pb + Pb collisions [23] . The solid curves show the SCM-based results while 
the dotted curve depicts the Gluon Saturation Approach [TT] . 
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